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Abstract. We study three-dimensional vortex lattice structures in purely dipolar 
Bose-Einstein condensate (BEC). By using the mean-field approximation, we obtain a 
stability diagram for the vortex states in purely dipolar BECs as a function of harmonic 
trap aspect ratio (A) and dipole-dipole interaction strength (D) under rotation. 
Rotating the condensate within the unstable region leads to collapse while in the stable 
region furnishes stable vortex lattices of dipolar BECs. We analyse stable vortex lattice 
structures by solving the three-dimensional time-dependent Gross-Pitaevskii equation 
in imaginary time. Further, the stability of vortex states is examined by evolution in 
real-time. We also investigate the distribution of vortices in a fully anisotropic trap 
by increasing eccentricity of the external trapping potential. We observe the breaking 
up of the condensate in two parts with an equal number of vortices on each when the 
trap is sufficiently weak, and the rotation frequency is high. 
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1. Introduction 

Vortices in Bose-Einstein condensates of alkali atoms first observed in the laboratory in 
1999 [[[]. Since then numerous experimental and theoretical studies along this direction 
were performed [2j. For instance, these studies include bending of vortex lines in a cigar 
shaped trap, array of orderly aligned lattices in the quantum-Hall regime, Tkachenko 
oscillations in the lowest Landau level and so on M- The rotational properties of 
BECs of alkali atoms are reviewed in Refs EE- Several techniques such as rotating 
the magnetic trap, laser stirring, decay of solitons, imprinting vortices using topological 
phases, superimposing an oscillating excitation to the trapping potential, and applying 
artificial magnetic fields mum were adopted in BEC experiments to nucleate the 
vortices. 

Early experiments and theoretical studies on vortices in BEC mostly focused 
on alkali Bose gases with local and isotropic interaction. In particular, there are 
many numerical studies proposed to analyse the stationary state of rotating alkali 
BECs [8]. However, the experimental realization of Bose-Einstein condensation in 
chromium ( 52 Cr) [HJ, dysprosium ( 164 Dy) jTUj . followed by erbium ( 168 Er) [TT] has 
enlightened new directions in understanding the properties of BEC in dipolar quantum 
gases. These recent works has revealed new phenomena due to the peculiar competition 
between isotropic short-range contact interaction and anisotropic long-range dipole- 
dipole interaction (DDI). The most significant features of the dipolar BECs are the 
emergence of biconcave shaped states, the dependence of the stability on the trap 
geometry, the roton-like dip in the dispersion relation for excitation waves, and the 
d- wave mode of the collapse P2UHEUIH]. The rotational properties of dipolar Bose 
gases have also been studied theoretically in the mean-field regime 

These studies revealed that the rotational properties of dipolar BECs are strongly 
influenced by the harmonic trap aspect ratio, DDI strength, contact interaction (Cl) 
strength and relative strengths between DDI and Cl [2D] EQ- Square and triangular 
vortex lattice structures have been predicted in dipolar BECs loaded in optical lattice 
potentials [22] . The stability of dipolar BEC is strongly dependent on the trap 
geometry and dipolar interaction strength. So far the regime for the stability has 
been analysed only for the non-vortex states of dipolar BECs without rotation [ 23 ;, 24 ]. 
Most of the previous studies on vortices in dipolar BECs are based on two-dimensional 
models only (THl HZ] HE] EH ED] • However, it will be more realistic to investigate the 
stationary vortex structures in full three-dimensions as they are readily comparable with 
experiments. 

In this paper, we study the vortex lattice structures of purely dipolar Bose-Einstein 
condensate by considering the full three-dimensional Gross-Pitaevskii equation. In 
particular, we analyse the stability regime for vortex state in purely dipolar BECs with 
respect to both trap aspect ratio and DDI strength. The stability of the vortices is 
confirmed by numerically evolving the vortex states in real-time. We show stationary 
vortex lattice structures for different trap aspect ratios within the stable regime. Further, 
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we study the collapse dynamics of biconcave shaped condensate during rotation. We 
calculate the chemical potential, rms (root-mean-square) radius and angular momentum 
of condensate as a function of rotation frequency. We estimate the number of vortices 
using Thomas-Fermi approximation and compare them with that obtained through 
numerical simulations. Finally, we notice the breaking of the condensate when the 
rotating dipolar BEC is in the fully anisotropic trap. 

The paper is organized as follows. In Sec. [2| we present the mean-field equation 
to study a dipolar BEC in a rotating trapping potential. In Sec. [3j we show the 
stability diagram for states without and with vortices in purely dipolar BECs as a 
function of harmonic trap aspect ratio and DDI strength. In Sec. [4j we study the 
stationary vortex structures observed numerically as well Thomas-Fermi calculation on 
some physical parameters during rotation of the condensate. Then we investigate the 
spatial distribution of vortices in a fully anisotropic trap in Sec. [5| Finally, in Sec. [6j 
we present a brief summary and conclusions. 


2. The mean-field Gross-Pitaevskii equation in rotating frame 


The pattern of vortices in the BEC can be studied using the mean field Gross-Pitaevskii 
(GP) equation [2| [25] . At ultra-low temperatures the properties of a dipolar Bose- 
Einstein condensate of N atoms, can be described by the mean-field GP equation in 
rotating frame with the nonlocal nonlinearity of the form [T6l 19]: 


i h 
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dt 
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where f dr|0(r, t) | 2 = 1. The trapping potential, V trap is assumed to be fully asymmetric 
of the form 

Vtrap(r) = ^m (calx 2 + uj 2 y 2 + u 2 z z 2 ) 

where a J x ,u) y and uj z are the trap frequencies, rri is the atom mass, and a is the atomic 
scattering length. L z = —i h(xd y — yd x ) corresponds to the 2 -component of the angular 
momentum due to the rotation of the dipolar BEC about 2 -axis with rotation frequency 
Q. 

For magnetic dipoles, the dipolar interaction is given by |26, 27] 

TT ro\ /h)/i 2 l-3cos 2 d /3cosV-l\ 

Udd{n) = 1V |R|3 (-2- ) ' (2) 

where R = r — r' determines the relative position of dipoles, 9 is the angle between R 
and the direction of polarization, /i 0 is the permeability of free space, and jl is the dipole 
moment of the atom, ip is the angle between the orientation of dipoles and 2 -axis. We 
consider the polarization of magnetic dipoles along the direction of 2 -axis as long as </? 
= 0 . 
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To compare the contact and dipolar interactions, it is convenient to introduce the 
length scale add = Pvh ‘ 2 ' m /(1 2 tt/ r 2 ) [28]. 52 Cr has a magnetic dipole moment of /7, = 6/ig 
(Hb is the Bohr magneton) so that add — 16ao, where ao is the Bohr radius. The 
dipole-dipole interaction strength is expressed as D — 3Na dd . 

Convenient dimensionless parameters can be defined in terms of a reference 
frequency Co and the corresponding oscillator length l = \JYij (mCo). Using dimensionless 
variables r ; = r/l, R/ = R /l, a' = a/l , a dd = a dd /l, t' = too, x' = x/l, y' = y/l , z' = z/l, 
= Q/ca, (f)' = 1'^ 2 4 k equation (|T|) can be rewritten (after dropping the primes from all 
the variables) as 


■ d<f>(r,t) 
1 dt 


— -V 2 + - (y 2 x 2 + u 2 y 2 + A 2 z 2 ) + 4TraN\(f)\ z 

Al £ 


-TtL z + 3Na dd I V d 3 d D (R)|0(r',t)| 2 dr' 
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with 

3 cos 2 <p — 1 A 
2 ) ’ 

r y = uj x /cu,u = co y /Cd,X = uj z /cu. We consider the cylindrically symmetric harmonic trap 
with 7 = v with u> x = co y = uo p and we use the reference frequency Co as cn p . For the 
fully anisotropic trap, the reference frequency is taken as the geometric mean, that is, 
Co = {(jj x 00y00 z ) C 3 , From now, we only refer to the dimensionless variables. 

We perform the numerical simulations of the 3D GP equation {(3j) using the split-step 
Crank-Nicolson method described in Refs. [22, 3RI. The dipolar integral in equation (|3]) 
diverges at a short distance in coordinate space. However, this can be circumvented by 
evaluating the integral in momentum space [23 EH- For the parameters we used in this 
work for the dipolar GP with (a = 0), there is no minimizer [32, 33], which means there 
is no actual ground state. So the states we refer here are in fact lowest local minimum 
states. 

The numerical simulations for the cylindrically symmetric and fully anisotropic 
traps are carried out with 192 x 192 x 192 and 320 x 128 x 128 grid sizes respectively, 
with Ax = Ay = A z = 0.2 (space step) and At = 0.004 (time step). 


V 3D 

bid 


(R) = 


1 — 3 cos 2 6 
(RU 


3. Non-vortex and Vortex states Stability 

The important feature of dipolar BEC is the emergence of unusual states. The lowest 
local minimum state looks like a biconcave shaped structure where the density maximum 
is not in the center of the condensate [231 El- Furthermore, the density oscillations of 
dipolar BEC with two and four peaks were also reviewed in a fully anisotropic trap. 
In this section, we are interested in studying the stability and structures in purely 
dipolar BEC in pancake trap. We have composed the stability diagram for a purely 
dipolar BEC as a function of the trap aspect ratio (A) and dipolar interaction strength 
(D). A similar stability diagram is given in Ref. [23]. A dipolar BEC is unstable and 
collapses for the number of atoms N above a critical value. This can be studied by 
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solving the three-dimensional time-dependent GP equation ([3j). We have broadened 
the stability diagram to A = 30 and calculated the stability region of the solutions. 
The biconcave shaped condensate density as shown in figure [IJa) is obtained for the 
parameters A = 7 and D = 30.4 from imaginary time propagation. In figure [l](b) we 
have shown the real-time evolution of the biconcave shaped state. One may note that 
similar biconcave structures are shown in Refs. Isa m- Also, we observed such local 



Figure 1. (Color online) Three-dimensional contour plot showing the transparent 
view of density oscillation of biconcave shaped condensate for A = 7 and D = 30.4 in 
the absence of rotation (f l = 0) evolve in real-time at time, (a) t = 0 and (b) t = 500. 
The contours levels are taken as \<j>(x,y, z)\ 2 = 0.03 and 0.04. 

density fluctuation for the following trap aspect ratios A = 7, 8,11,12,15,16,19, and 
20. The structured condensates are locally stable in the sense that they are stable only 
within a local minimum of the energy. Further increase in the number of particles or 
increase in dipolar interaction energy will change the state to unstable. Here we are 
interested in rotating this biconcave shaped condensate to study the time evolution. 
With this objective, we prepared the initial state solution by solving the equation 0 
using imaginary time propagation in the absence of rotation (Q = 0) and then seed 
the resulting solution in real-time propagation by applying a rotation with frequency 
Q = 0.5. At the time of rotation, the dipoles are immensely pulled towards the outer 
rim, and two peaks formed as shown in figure [2|R). These two peaks sustain for a very 
short time, and biconcave shaped condensate becomes unstable. 

If we analytically evolve an exact solution with zero angular momentum in the z 
direction, then term —QL Z in equation 0 plays no role. However, this may not be 
the case either physically or if it is evolved numerically. Physically, there is always 
some fluctuation in the wave function field. Fluctuations are random and have in 
principle contributions of all components of different angular momentum. In the 
numerical calculations, roundoff errors are always present and give nonzero angular 
momentum contribution. The term —QL Z will then act on the perturbations, growing 
them exponentially (modulational instability) and the condensate eventually collapses. 
The collapse dynamics during rotation is shown as two- and three- dimensional view in 
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Figure 2. (Color online) The condensate density y, 0)| 2 for A = 7 and D = 30.4 
in xy plane showing the collapse dynamics of rotating biconcave shaped condensate 
with Q = 0.5 in real-time at time, (a) t = 0, (b) t = 60, (c) t = 64, and (d) t = 65. 


figures [2] and [3j respectively. Further, we have constructed a phase diagram illustrating 



Figure 3. (Color online) Three-dimensional contour plots of the density |0(x, y, z)\ 2 
for A = 7 and D = 30.4 showing the rotating biconcave shaped condensate with ft = 0.5 
in real-time propagation at times (a) t = 0, (b) t = 60, and (c) t = 64. The contour 
levels are taken as \<j)(x, y, z)\ 2 = 0.04. 

the stability region for stable vortex lattice in A —D parameter space as shown in figure [4} 
The stable region for vortex lattice is located below the stability region of the states 
(dashed blue line with empty circles) in the figure [4j To carry out this observation, we 
have prepared the solution using imaginary time propagation with Q = 0 and progress 
the solution in imaginary-time with rotation frequency Q = 0.7. Nevertheless, the 
stability diagram has been checked for the rotation frequencies in the range 0.7 and 0.99 
and no significant changes in the stability boundary are observed. Also, we examined the 
stability of stationary vortex lattice solutions by evolving in real-time. In the figure [4| 
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Figure 4. (Color online) Stability diagram of a purely dipolar BEC in cylindrically 
symmetric trap as a function of trap aspect ratio (A) and dipolar interaction parameter 
(. D ). Below the dotted blue line the local minimum with non-vortex state is stable while 
below the solid red line is the region of stable vortex lattice. 


the region below the solid red line with empty circles corresponds to stable vortex 
lattice. The vortex states are not stable in the region between the dashed blue with 
empty circles and solid red line with empty circles. In the next section, we discuss the 
feasible vortex structures in the shaded region below the solid red line. For the present 
study, we mainly consider two different harmonic trap aspect ratios, A = 10 and 30, and 
the dipolar interaction strength is chosen within the vortex lattice region in the figure [4| 

4. Stationary Vortex lattices in purely dipolar BECs 

In the following, we show several stationary vortex structures in purely dipolar BEC in 
the different harmonic trap aspect ratios A = 10 and 30, and correspondingly the dipolar 
interaction strength is chosen as D ~ 38 and 300, respectively. First, in figure [5j we 
show the stationary vortex structures in A = 10. We prepare the initial wave function 
by solving equation (§ in the absence of rotation (Q = 0). The vortices are then created 
by using imaginary time propagation with the inclusion of rotation (12 ^ 0) to observe 
the stationary vortex structures. When the condensate begins to rotate the multiply 
quantized vortices enter into the condensate from the surface. As time progress, these 
vortices approach to a stationary vortex configuration as shown in figure [5} We notice 
the basic configuration of centered single vortex surrounded by five vortices at rotation 
frequency 12 = 0.7 as shown in figure [5](b) . At 12 = 0.8 the vortex structure is shown 
in figure [5](c) resembles the structures observed in the references [16, 20| for dipolar 
condensates. Further, we observed the square lattice at 12 = 0.9. To test the stability 
of the vortex structures, the relaxed solution is evolved in real -time for the same 
parameters of imaginary time. The vortex structure persisted as depicted in figure [6j 
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Figure 5. (Color online) Three-dimensional view of condensate density \c/)(x,y, z)\ 2 
with vortex lattice of purely dipolar BECs with add = 16 ao, a = 0, D = 38, A = 10 
at rotation frequencies (a) Q = 0.5, (b) Q = 0.7, (c) Q = 0.8 and (d) Q = 0.97. The 
contour levels are taken as |0(x, y, z)\ 2 = 0.0125 
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Figure 6. (Color online) Two-dimensional view of the stable evolution of the 
condensate density |0(x, ?/, 0)| 2 with vortices arranged in a square lattice with add = 
16ao, a = 0, A = 10, D = 38, and Q = 0.9 at time, (a) t m 0 and (b) t — 800. 


even after long time evolution. 

At high rotation frequency (e.g. El — 0.97), we found 16 vortices as shown in 
figure [5](d) . Moreover, one may note from figure |4] that, for obtaining the vortex lattice 
in the condensates with larger dipolar strength, it is necessary to consider significantly 
larger A. In this case we choose the trap aspect ratio A = 30 with dipolar interaction 
strength D ~ 300 corresponding to about 118 000 52 Cr atoms. We observed a similar 
arrangement in the vortex structures in A = 10 and 30, when the number of vortices are 
equal. It tells ns that the condensate tries to persist with the similar vortex structures, 
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Figure 7. (Color online)Three-dimensional view of condensate density \4>(x, g, z)\ 2 
with vortex lattice of purely dipolar BECs with add = 16 ao, a = 0, D = 300, A = 30 
at rotation frequencies (a) f2 = 0.27, |0(x, y, z)\ 2 = 0.0125, (b) Q = 0.3, |</>(x, y, z)\ 2 = 
0.005, (c) 12 = 0.4, \(j)(x,y, z)\ 2 = 0.005 and (d) 12 = 0.5, \c/)(x , y , z)\ 2 = 0.005 (e) 
12 = 0.8, \cj)(x,y,z)\ 2 = 0.005, and (f) 12 = 0.97, |0(x, y, z)\ 2 = 0.004. 

even though the trap aspect ratio or dipolar interaction strength are distinct. For 
instance, we noticed a qualitatively similar array of vortex lattices for both A = 10 (see 
figure [6]) and A = 30 (see figure [7|(b)) , with fi = 0.9 and 0.3, respectively. We observe 
qualitatively similar structures for A = 10 (see figure [5](d) ) and A = 30 (see figure gc)), 
with C = 0.97 and 0.4, respectively. The condensate in a strong pancake trap creates a 
large number of vortices, even at low rotation frequency. 

Further, we observed some distortion in the vortex lattice structure when the 
number of vortices is sufficiently large. In particular, the distortion is larger near to 
the surface than near to the center of the condensate. All the vortices have the same 
charge of vorticity, and their repulsion keeps the stable configuration. We note that 
the size of the vortex core radius becomes significantly larger in the vicinity of the 
surface than in the middle of the condensate while increasing the rotation frequency. 
Consequently, near to the surface, repulsive interaction diminishes the density of vortices 
and leads to the distortion in the vortex lattice. It is equivalent to observing transverse 
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shear waves in quantum gases for large rotation frequencies jl]. In this regime, it takes 
place a decrease in the elastic shear strength of the vortex lattice. 

In figures [8] and [9j we plot the calculated values of chemical potential, rms radius, 
and angular momentum as a function of rotation frequency for A = 10 and 30. First, 
we calculate the chemical potential with respect to the rotation frequency in mean- 
held Thomas-Fermi (TF) regime. When the interaction energy is large compared to 
the kinetic energy, the kinetic energy can be neglected and enters into TF regime. We 
assume the normalized density of the dipolar BEC of the form [3H 35, 36, 137] 


n(r,f) = |</>(r, f)| 2 = 


15 
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v2 1 


8*mmM l m) 
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where R p (t) and R z (t) are the radial and axial sizes. In the TF regime one has the 
following set of coupled ordinary differential equations for the evolution of the condensate 
sizes [35]: 
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Here the atomic scattering length is taken as a = 0 so that the dipolar interaction shows 
the dominant effect. It may be noted that the external rotation expands the condensate 
radially and shrinks it axially. As a consequence, the dependence of TF radii on Q can 
be given by [2j, 


R P ( fi) 
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R P { o) v 7 ’ RM 

Also, the chemical potential has the form, 

= Ptf{ 0) (l — fl 2 ) ^ • 


= (i - « 2 ) 
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( 10 ) 


As shown in figure |8|a) the chemical potential decreases continuously and it goes 
to zero when Q = 1.0. In figure 8[b-e) we plot the rms radius, \J (r 2 ) = 

(J cf)*(r,t) r 2 b(r. t) dr) , and chemical potential as a function of rotation frequency 
calculated from the numerical solution of the GP equation. As expected, the condensate 
rms radius increases with the increase of rotation frequency due to the expansion of 
condensate. The rms radius is shown in figures [8|(b) and (c) for A = 10 and 30, 
respectively. One may note that, in the absence of rotation, the numerically calculated 
chemical potential compares well with TF results [31]. Whereas, in the presence 
of rotation, the TF chemical potential is about two times less than the numerically 
calculated value. In figures |8](d) and (e) we show the variation of chemical potential for 
A = 10 and 30 respectively. Next, we calculate expectation value of angular momentum 
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Figure 8. Plot of the chemical potential and rms radius as a function of the rotation 
frequency Q for the dipolar BEC with the scattering length a = 0, A = 10 with D ~ 38 
and A = 30 with D ~ 300. (a) TF chemical potential calculated from Eqn. ( [To] ) . Plot 
of the numerical calculation of radius and chemical potential in (b), (c) and (d), (e) 
respectively. 



Figure 9. (color online) Plot of the expectation value of angular momentum ( L z ) as 
a function of the rotation frequency for the purely dipolar BEC trapped in A = 10, 
30 and 38, 300, respectively. 
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as a function of 12 from the numerical solution. Figure [9j depicts the plot of (L z ) with 
respect to 12 for A = 10 and 30. The increase in the (L z ) is associated with the entry 
of vortices into the condensate. We are interested in studying the dependence of the 
vortex number (N v ) on the rotational frequency 12 in purely dipolar BEC. The rotating 
condensate has a dense array of vortices with a uniform density n v = rn{iluj)/Till. The 
number of vortices present in the condensate is [ 2 ], 


N v = 


miThsj) 

h 


R 2 M- 


(ii) 


The number of vortices increases linearly with 12 assuming the BEC in a fixed axially 
symmetric harmonic trap. The radius R p { 0) in TF regime can be calculated from the 
coupled equations ([ 6 ]) and (| 7 |), and the corresponding i? p ( 12 ) can be obtained from 
equation ([9j) . In figure [To|(a) , we have estimated the number of vortices in the different 
harmonic trap aspect ratios A = 10, 20 and 30 in TF regime and compared with 
numerically calculated equilibrium numbers. We noted some deviations in the number 
of vortices from TF results and calculated numerically. To compare the deviation in the 
number of vortices in BEC without DDI, we calculate N v for the BEC without DDI by 
tuning to the magic angle tp = 54.7°, where the dipolar interaction averages to zero in 
equation Q. Then we obtain N v at a = 50ao, 100ao and we observe the deviation as 
shown in figure [Toj^b) . In particular, when 12 > 0.9 the deviation is larger and it is due 



11 12 

Figure 10. (color online) Plot of the equilibrium vortex number (N v ) as a function of 
the rotation frequency fl for the (a) purely dipolar BEC trapped in A = 10, 20 and 30, 
with D = 38, 160 and 300, respectively and (b) non-dipolar BEC trapped in A = 30. 

to the strong expansion of condensate during the rapid rotation. 

5. Vortex lattice in fully anisotropic trap 

Finally, we concerned the purely dipolar BEC in fully anisotropic trap with 7 ^ v 7 ^ A 
in the equation (|3j). Whenever 7 //, the system breaks the cylindrical symmetry. For 

the present analysis, we fix A = 10, v = 1, and vary 7 to study the changes in the 
spatial distribution of vortices in purely dipolar BECs by increasing the eccentricity 
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Figure 11. (Color online) Vortex lattices in rotating dipolar BECs trapped in 
asymmetric trap with 7 = 0.5, A = 10, D = 38 and rotation frequency: (a) Q = 0.425, 
(b) Q = 0.6, (c) £2 = 0.7. The contour levels are taken at |0(x, y , z) | 2 = 0.01, 0.005 and 
0.002, respectively. The corresponding phase of the condensate density |0(x,t/,O)| 2 is 
shown in d-f. Blue lines are density contour lines. Blue circles show the position of 
vortices, red squares show the center of hidden vortices. 

of the trapping potential along x-direction. The condensate strongly elongates along 
the x-direction at 7 < 0.5. The centrifugal force also elongates the condensate along 
x-direction. We have shown four linearly arranged vortices in figure [7T|(a) for 7 = 0.5 
and £2 — 0.425. In the figures 0 and [ 12 J we observed linear and zig-zag arrangements 
of vortices with respect to the strength of eccentricity of the trap, respectively. Such 
linear and zig-zag vortex configurations have been observed in conventional quasi two- 
dimensional BEC [40J • Further increase in Q splits the condensate at a critical rotation 
frequency (Q sp ). In addition to the repulsive dipolar force between condensate atoms, 
the rotation provides supplementary elongation at fl > fl sp . Splitting is a combined 
effect of rotation and dipolar forces. As rotation increases, the condensate creates extra 
vortices and also stretches. In the limiting case, it remains a vortex in the center 
and further stretching breaks the condensate in two parts. With further increase of 
rotation, no more vortices are generated since extra angular momentum can now be 
accommodated by orbital angular momentum. Stretching is enhanced depending on 
the trap geometries. We have observed only two vortices as shown in figure [l2|a) at 
7 = 0.7 for the rotation frequency fl = 0.425. This is exactly half the number of vortices 
when compared to that with 7 = 0.5 for the same rotation frequency. This evidences 
the dependence of the number of vortices on the eccentricity of the trap. At sufficiently 
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(a) (b) (c) 



Figure 12. (Color online) Vortex lattices in rotating dipolar BECs trapped in 
asymmetric trap with 7 = 0.7, A = 10, D = 38 and rotation frequency: (a) Q = 0.425, 
(b) = 0.7 , (c) Q = 0.91. The contour levels are taken at |0(x, y, z)\ 2 = 0.02, 0.012 

and 0.002, respectively. The corresponding phase of the condensate density |0(x, y , 0 )| 2 
is shown in (d-f). Blue lines are density contour lines. Blue circles show the position 
of vortices, red squares show the center of hidden vortices. 


high the condensate splits into two parts as shown in figure [l2](c) . Here the splitting 
occurs at Q sp = 0.91, which indicates that the trap with large 7 requires substantially 
larger rotation frequency than the condensate with smaller 7 . The phase of the vortices, 
marked with circles, is shown in figures ©d)- (f) and [T 2 p)-(f). We also marked few 
more vortices as red squares. These vortices are called hidden (or ghost) vortices since 
they cannot be observed in 3 D density distribution. Due to hidden vortices, there is 
a phase defect distributed along the y -axis in the split condensate. A similar phase 
distribution is observed in rotating BEC in a double-well potential [3.8]. The splitting 
with vortices are also observed in normal BEC, where it is triggered by random noise or 
a complex frequency to the stationary vortices [39]. However, splitting is spontaneous 
in rotating dipolar BEC. 

we show fl sp at different trap aspect ratios, A = 10, 30, 50 and 


In figure 13 


100 with respect to the eccentricity of the trap 7 . When the anisotropy is fixed as 
7 = 0.5 and A = 10 then the splitting is observed at fl sp = 0.65. One can observe that 
Cl sp increases with increasing trap aspect ratio. For instance, when A = 50 and 100 
the corresponding required fl sp are 0.72 and 0.78, respectively. The strong axial trap 
preserves the condensate from splitting, and one needs to apply stronger fl to observe 
it. On the other hand, if we increase the anisotropy to 7 = 0.7, the Q sp also increases 
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Figure 13. Critical rotation frequency (fl 3p ) for the splitting of the condensate 
with respect to anisotropy ( 7 ) of the trap along ^-direction for the parameters v = 1 , 
D = 38 and for different trap aspect ratios A = 10, 30, 50 and 100. 


to 0.9 for A = 10. The condensate does not split in strong enough axial traps as for 
7 = 0.7 and A = 50 and 100 . We have also carried out the calculations with an effective 
2 D GP equation in a rotating frame and observed a similar splitting with same vortex 
distribution for the corresponding parameters used in 3D calculations, confirming our 
results. 

6. Summary and Conclusion 

In summary, we have studied the vortex lattice structures in purely dipolar Bose-Einstein 
condensate of 52 Cr atoms by considering the full three-dimensional Gross-Pitaevskii 
equation. We have identified the stability regimes for non-vortex as well as vortex 
states in purely dipolar BECs with respect to both trap aspect ratio and dipole-dipole 
interaction strength. The stability of vortex lattice structures are confirmed by real 
time evolution. We have shown stationary vortex lattice structures for different trap 
aspect ratios within the stability regime. Also, distortions in the lattices at high rotation 
frequencies have been observed. 

We estimated the number of vortices using Thomas-Fermi approximation and 
compared with those from numerical simulations. Further, we have analysed the spatial 
distribution of vortices in the fully anisotropic trap and observed linear and zig-zag 
arrangement of vortices. We noticed the breaking up of the condensates into two 
parts with an equal number of vortices. The splitting occurs due to a combination of 
centrifugal force due to rotation and repulsive dipolar interaction. Increasing rotation 
frequency creates extra vortices and stretches. Close to tt sp it appears a vortex at 
the center and further stretching breaks the condensate in two parts. Critical rotation 
frequency for splitting depends on the trap aspect ratio and dipolar parameters. 
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The predicted phase diagram of the stable vortex state will be useful to demonstrate 
the parameters such as dipolar strengths and trap aspect ratios for making experimental 
and theoretical studies on rotating dipolar quantum gases. Further, the stability regime 
is relevant to investigate rapidly rotating dipolar BECs in the lowest Landau level. We 
found that dipolar BECs under rotation can produce square vortex lattices. This was 
found with two-component BECs m and recently in dipolar Fermi gas 021 • Observing 
the breaking in the rotating dipolar BECs in a fully anisotropic trap will be a new 
experimental exploration. 
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